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Abstract 

This paper is part of a study of the frictional dynamics of a confined solid lubricant 
film - modelled as a one-dimensional chain of interacting particles confined be- 
tween two ideally incommensurate substrates, one of which is driven relative to the 
other through an attached spring moving at constant velocity. This model system 
is characterized by three inherent length scales; depending on the precise choice of 
incommensurability among them it displays a strikingly different tribological behav- 
ior. Contrary to two length-scale systems such as the standard Frenkel-Kontorova 
(FK) model, for large chain stiffness one finds that here the most favorable (lowest 
friction) sliding regime is achieved by chain-substrate incommensurabilities belong- 
ing to the class of non-quadratic irrational numbers (e.g., the spiral mean). The 
well-known golden mean (quadratic) incommensurability which slides best in the 
standard FK model shows instead higher kinetic-friction values. The underlying 
reason lies in the pinning properties of the lattice of solitons formed by the chain 
with the substrate having the closest periodicity, with the other slider. 
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1 Introduction 



Nonlinear systems driven far from equilibrium exhibit a very rich variety of 
complex spatial and temporal behavior. In particular, in the emerging field of 
nanoscale science and technology, understanding the nonequilibrium dynamics 
of systems with many degrees of freedom which are pinned in some external 
potential, as commonly occurs in solid-state physics, is becoming more and 
more a central issue. Classic textbook examples are the motion of disloca- 
tions in metals, of domain walls in ferroelectrics, the problem of crowdion in a 
metal, sub-monolayer films of atoms on crystal surfaces, transport in Joseph- 
son junctions, the sliding of charge density waves and so on. Sliding friction of 
solids belongs to this class of systems as well, because the microscopic periodic 
asperities of the mating surfaces may interlock. 

One of the pervasive concepts of modern tribology with a wide area of relevant 
practical applications as well as fundamental theoretical issues is the idea of 
free sliding connected with incommensurability. When two crystalline work- 
pieces with lattices that are incommensurate (or commensurate but not per- 
fectly aligned) are brought into contact, the minimal force required to achieve 
sliding (i.e., the static friction, F s ) should vanish, provided the two substrates 
are stiff enough. In this configuration, the lattice mismatch causes the total 
energy to be completely independent of the relative position of the sliders, 
whereas hardness can prevent pinning and the associated stick-slip motion of 
the interface atoms, with a consequent negligibly small frictional force. The 
remarkable conclusion of frictionless sliding can be drawn, in particular, in 
the well explored ID context of the Frenkel-Kontorova (FK) model (see [1] 
and references therein). Experimental observation of this kind of superlubric 
regime has recently been reported [2]. 

Based on the success of this kind of simplified tribological approach [3], the 
hope is to be able to learn more from the direct study of microscopic dy- 
namics in similarly simplified ID models. The substrates defining the sliding 
interface are modelled as purely rigid surfaces or as ID (or 2D) arrays of parti- 
cles interacting through simple (e.g., harmonic) potentials. Despite its extreme 
idealization, this kind approach seems able to anticipate features of the actual 
experimental results or of more complex molecular dynamics simulations of 
frictional phenomena. In this context, the application of driven FK like models 
(see [4] and references therein), describing the dissipative dynamics of a chain 
of interacting particles that slide over a rigid periodic substrate potential due 
to application of an external driving force, has found an increasing interest as 
a possible interpretative key to understand the atomic processes occurring at 
the interface of two materials in relative motion. The essential feature of the 
static and dynamic properties of the FK model consists in the competition 
between the interparticle interaction (having natural lattice constant b) and 
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the substrate periodic potential (of spatial period a). Indeed, if the former 
favors a uniform separation b between particles, the latter tends to pin the 
atom positions to the bottom of the wells, evenly spaced by the period a. This 
length-scale competition, often referred to as frustration, results in a fasci- 
nating complexity of spatially modulated structures. For this two competing 
length-scale system, following the pioneering work of Aubry [5], the tribolog- 
ical properties depend mainly on the value of the winding number a/b being 
rational (commensurate interface) or irrational (incommensurate interface). 
Rational interfaces are always ingrained and frictional, whereas irrational in- 
terfaces exhibit an "Aubry transition" from frictional to frictionless as the 
chain stiffness increases. 

In real situations, however, such a case of "dry" , unlubricated friction is often 
exceptional. The physical contact between two solids is generally mediated 
by so-called "third bodies", which may act like a lubricant film. If the third 
body is a solid lubricant, and in particular a crystalline one, the sliding inter- 
face corresponds in fact to a system with three inherent lengths: the periods 
of the bottom and top substrates, and the period of the lubricant layer. As 
sketched in Fig. 1, in order to study the role of incommensurability among 
the three interface inherent lengths on the sliding dynamics, we may consider 
a ID generalized FK model consisting of two rigid sinusoidal substrates, of 
spatial periodicity a (bottom) and c (top), harassing a chain of harmonically 
interacting particles, of equilibrium length b, mimicking the confined lubri- 
cant. Guided by a similarity with the recent study of the driven quasiperiodic 
FK model [6], we show that, for sufficiently stiff chain (hard crystalline lu- 
bricant), the best low-friction regime is achieved for incommensurabilities a/b 
and c/a related through non-quadratic irrational numbers (as, for example, 
the cubic spiral mean [7]) rather than through quadratic irrationals (e.g., the 
golden mean). This result demonstrates how, in a t/iree-length sliding sys- 
tem, the kind of interface incommensurability (and not just the distinction 
between commensurate and incommensurate geometry) can dramatically in- 
fluence, both quantitatively and qualitatively, the tribological behavior. 



2 Model and kinetic friction 

In numerical simulations of frictional processes, different approaches may be 
used to induce motion: a constant- velocity algorithm, where the top substrate 
is directly forced to slide with a given velocity, or a constant-force algorithm, 
where the external driving force is applied to the top substrate (or even directly 
to the lubricant layer), or finally a spring-force algorithm, where the slider is 
driven through an attached spring moving at constant speed. The latter pulling 
procedure can be viewed as a way to mimic not only the experimental driving 
device, but to some degree also the elasticity of the slider. 
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In this paper, we consider a ID system of two rigid sinusoidal substrates and a 
chain of interacting particles embedded between them [8], as shown in Fig. 1. 
The top substrate of mass M is driven through a spring, K ext , which is pulled 
with constant velocity V cxt . The equations of motion become 
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where x^ (i = 1, . . . , N) and X top stand for the coordinates of the N chain 
particles and the top substrate, respectively. The damping 7-terms describe 
the dissipative forces which are proportional to the relative velocities of the 
lubricant atoms with respect to both rigid substrates. The viscous coefficient 7 
represents degrees of freedom inherent in the real physical system (such as sub- 
strate phonons, electronic excitations, etc.) which are not explicitly included 
in the model. The last (sinusoidal) terms represent the on-site interaction be- 
tween the particles and the substrates. The amplitudes of these rigid potentials 
are chosen such that the same factor (1/2) sits in front of their derivatives in 
the equations of motion - whence the periodic forces they produce have the 
same magnitude. The interparticle chain interaction [fourth term in Eq. (1)] 
is harmonic with strength K and equilibrium spacing b. We use dimensionless 
units with chain atom mass m — 1 and a bottom substrate period a — 1. 
We focus here on the tribologically most interesting case of the underdamped 
regime, where the damping coefficient 7 is much smaller than both character- 
istic vibrational frequencies of one particle sitting at the minima of the two 
substrate potentials. 



The system is initialized with the lubricant particles placed at rest at uniform 
separation b. After relaxing the starting configuration, the stage attached to 
the top substrate through the spring K ext , begins to move at constant speed 
V cxt . The equations of motion (1), (2) are integrated using the standard fourth- 
order Runge-Kutta algorithm. After reaching the steady state, the relevant 
physical characteristics of the system are measured. The detailed behavior 
of the driven system in Eqs. (1), (2) depends crucially on the length-scale 
relative ratios a/b and c/a of the substrates and the chain. Since commen- 
surability can only be achieved under well-controlled experimental operating 
conditions [9], in this investigation we focus on a few specific but typical cases 
of incommensurability among the three length scales a, b, c. Specifically we as- 
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sume a single ratio a/b = c/a and consider the two irrational cases previously 
studied in the context of the driven quasiperiodic FK model [6], namely the 
golden mean and the spiral mean [7]. Nevertheless, it should be noted that the 
qualitative features of the observed dynamics (in a similar ID confined model 
corresponding to the M — > oo limit of Eq. (2), so that in Eq. (1) X top = V^t) 
have been shown recently to survive for much more general values of a/b and 
c/a [10,11,12], so that this specific choice of incommensurability should not 
be considered too restrictive. By imposing periodic boundary conditions in 
order to simulate an infinite system, we are forced to approximate the desired 
irrational winding numbers a/b, c/a by ratios of integers. Continued-fraction 
expansion [13] provides a standard and well-established technique to serve this 
purpose. Quantitatively, a given rational approximant can be checked to be of 
sufficiently high order to model the desired incommensurability. The chain of 
iV atoms and length L = Nb is thus confined in between the bottom and top 
substrates with N a = L/a and N c = L/c minima, respectively. 



In tribology applications, one of the main issues is the kinetic friction force F^ in 
which describes the energy losses inside the contact layer when the substrates 
move relatively to each other. In our model, F^ n can be easily evaluated from 
energy balance. When the top substrate moves for a distance AX top = X top At 
relative to the fixed bottom substrate, due to the friction force the system 
dissipates an energy 
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In the regime of steady motion, this must balance the losses due to the dissi- 
pation generated by the 7-damping terms in the chain, i.e., 
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where At = t 2 — t±. Thus, for the kinetic friction force, we obtain 
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In order to highlight in the kinetic friction force the dissipative part resulting 
from the internal dynamics of the confined chain, we subtract the systematic 
contribution arising from the motion of the top substrate, and also normalize 
Fkin to the number of particles N, i.e., 

/kin = ~ ^T^top- (6) 
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3 Results 



In the range of model parameters considered, our numerical simulations reveal 
a strikingly different sliding behavior of the system depending on incommen- 
surability. In particular, it is not sufficient simply to distinguish between com- 
mensurate and incommensurate periodicities. The "degree" of incommensura- 
bility (as measured by continued-fraction expansion) plays a fundamental role 
as well. Here the mathematical properties of the length-scale ratio leave the 
realm of abstract number theory to become physically relevant. In the context 
of the M — > oo confined model, the peculiar incommensurability-dependent 
phenomenology has been recently explained [10] by the top slider (meaning 
here that one among the two substrates whose period c is the longest, so that 
b < a < c) rigidly "dragging" the topological solitons, with linear density 
p s — 1/6 — 1/a, that the chain forms with the bottom substrate. 

In particular, the choice of incommensurability among the three length scales 
a, b, c governs the sliding behavior by inducing a commensurate or incommen- 
surate matching of the soliton lattice in the chain with respect to the top 
substrate. In the golden-mean case, p s c = 4> 2 — <fi = 1, the soliton lattice is 
commensurate with the top slider, and so it must ingrain with it. In the spiral 
mean case, p s c = a 2 — a is irrational, the soliton lattice is incommensurate 
with the top slider, and so it can slide freely over it. That describes precisely 
what happens [14]. 

Figure 2 shows the kinetic friction force /kin as a function of the interparticle 
interaction strength K for different values of the velocity of the moving stage. 
The plots are given for two irrational choices of the three length scales: the 
quadratic golden mean (N a = 144, N = 233, N c = 89) and the cubic spiral 
mean (N a = 265, N = 351, iV c = 200). By simulating larger systems, we 
checked that these rational approximants [7] are of sufficiently high order to 
mimic the selected incommensurabilities. In all simulations we set M = N c m, 
and take 7 = 0.2. 

For the golden incommensurability, the qualitative (monotonic) behavior of 
the kinetic friction does not depend significantly on either K or V cxt . For the 
same value of the chain stiffness, /km decreases with decreasing driving velocity. 
The chain dynamics exhibits an asymmetric sliding with respect to the two 
sinusoidal substrates, moving as expected with an intermediate mean velocity, 
but always faster with respect to the top substrate with the longer spatial 
period. We checked that the asymmetry persists even when the amplitudes 
of the substrate potentials are varied significantly with respect to each other. 
This robust sliding state characterized by a regular time-periodic dynamics, as 
demonstrated in Fig. 3(b), is specific to the golden incommensurability and to 
all other quadratic irrationals that we considered. Each single particle in the 
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chain performs exactly the same periodic motion, differing simply by a phase 
shift. The period T m 26.180 of the motion is highlighted by the corresponding 
frequency peak v — (1 — -) V cxt ~ 0.038197 in the Fourier transform of the 
velocity Xi of one generic particle (Fig. 3(c)). 

As anticipated above, specific quadratic irrationals A, such as the golden mean, 
occur to satisfy p s c = A 2 — A = m/n, so that the soliton lattice and the top 
slider are commensurate. This commensurability provides an intrinsic ingrain- 
ing which extends to in principle arbitrarily large chain stiffness, without any 
transition. As an additional rationale, one may imagine that this particularly 
strong asymmetric center of mass motion is precisely such as to excite, via its 
time-periodic feature, parametric resonances inside the chain [15], thus con- 
verting part of the center-of-mass translational kinetic energy into internal 
vibrational excitations, making in turn the golden irrationality less favorable 
to sliding than the spiral one at the values of K considered here. 

For inherent lengths related instead by the cubic spiral mean, the function 
fkin(K) is found to display a non-monotonic and richer behavior, with sev- 
eral different regimes. Only at high driving velocities, V cxt > 0.01, these ex- 
otic features diminish and eventually disappear. The single-particle motion 
(see Fig. 3(a)) of the chain, as confirmed by the Fourier spectrum of Xi (not 
shown), is now definitely not periodic. This suppresses the parametric reso- 
nances inside the confined layer as possible channels for energy dissipation, 
making, for stiff chains (large K), the spiral-ratio slider more effective than 
the golden-ratio one. As noted above, in the spiral case the soliton lattice is 
incommensurate with the upper slider, and when stiff enough it may actually 
slide essentially freely over it. The evolution at large K seems in fact towards 
an almost frictionless smooth sliding. For smaller stiffness, a sharp change 
in the sliding behavior of the spiral chain is observed at a critical stiffness 
K c ~ 5.2, near the value where the Aubry (static depinning) transition for the 
quasiperiodic FK model was found previously [6]. When the driving velocity 
decreases (see lower curves in Fig. 2), a plateau region of larger dynamic fric- 
tion f\a n {K) develops, in sharp contrast to the unpinning observed in the static 
behavior. As can be seen, the location of its left edge is almost independent 
of V cxt and centered around the value K c . For K < K c finally the dynamics 
is characterized by a stick-slip motion. This is consistent with the forced slid- 
ing of an incommensurate, but slack and therefore pinned, soliton irregular 
arrangement. At large and increasing chain stiffness, the spiral friction drops 
considerably (with a center of mass velocity becoming symmetric, i.e. V cxt /2), 
as compared to the golden-mean case: in the latter, friction retains quite large 
values due to the velocity asymmetry of its fully regular time-periodic motion 
which is not destabilized until strikingly large K. 

In order to get more insight into different regions to the left, in between, 
and to the right of the spiral frictional plateau, we analyze in detail the four 
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points circled in black in Fig. 2, all belonging to the same driving velocity 
Kxt — 10~ 4 . The four upper panels of Fig. 4 show the time evolution of the 
spring force F ext = K cxt (X top — V ext t); the lower panels display the correspond- 
ing velocities of the top substrate and of the center of mass of the chain. The 
first upper panel for K = 4.0 clearly displays a typical sawtooth dependence 
of F ext (t), the hallmark of stick-slip motion. Stick-slip as a phenomenon is 
quite general and extremely well known. It has been observed in many simu- 
lations of microscopic friction, as well as in the tip dynamics of atomic force 
microscopy and in surface force apparatus experiments of confined lubricants 
under shear (see, e.g., Ref. [9] and references therein). Energy is stored in 
the springs while atoms are trapped in a metastable state (sticking), and it 
is converted to kinetic energy as the atoms pop to the next metastable state 
(slipping). The value of the spatial distance covered during the slip events 
and the regular-periodic or intermittent-chaotic temporal alternation between 
stick and slip can be affected by several factors: the intrinsic periodicity of 
the lattices defining the substrate potentials, the potential energy landscape 
of accessible metastable states, the specific features of imposing driving, etc. 
At K = 4.6 the slip events are prevailing, but the stick-slip regime of the chain 
is still clearly visible. The next value examined, K = 6.6, lies in the plateau 
region. Here the chain stiffness is too high to support stick-slip, and the con- 
fined layer motion reaches a sliding regime with almost constant value of the 
pulling spring force, and the velocity of the top substrate exhibiting only very 
tiny oscillations around V ext - The increase in the average kinetic friction /km of 
the plateau states (relative to the contiguous stick slip dynamics region) can 
be ascribed to the asymmetric chain motion with respect to the two substrates 
in this regime. In this range of stiffness and driving velocity values, a better 
dynamical phase-matching is possible between the confined layer and the two 
sinusoidal substrates, promoting the excitation of dissipative parametric res- 
onances, similar to the golden case. Inside this plateau, even if the detailed 
particle dynamics still remains irregular (quasiperiodic), a remarkable phase 
cancellation between the Fourier spectra of different chain particles (all having 
the same amplitude spectrum, with different phases) yields strictly periodic 
center-of-mass motion of the embedded chain. Finally, at K — 8.6 the con- 
fined layer slides almost freely with a very low value of /kin- Sliding becomes 
symmetric in this regime, i.e. the average chain velocity nearly equals one half 
of the relative velocity of the two rigid substrates. 



4 Conclusions 



When two periodic substrates are forced to slide relative to each other but 
a solid "lubricant" layer is interposed in between, the frictional dynamics of 
the intermediate layer may display an unexpected behavior, significantly de- 
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pendent on the incommensurability chosen. Our results show that, depending 
on physical parameters such as the lubricant stiffness, at least for a 3-lengths 
contact interface (a typical physical situation) there is no quantitative and 
qualitative uniformity of behavior in incommensurate sliding friction, and that 
certain incommensurate (e.g., non-quadratic) geometries can favor sliding sys- 
tematically better than others. In particular, we showed here that for suffi- 
ciently stiff confined chains, the standard golden-mean case reveals a higher 
kinetic friction than the spiral one. This is attributable to the hidden commen- 
surability of the soliton lattice formed between the chain and the lower slider, 
which occurs in the golden mean but not in the spiral mean case. It is not 
clear at the moment to what extent this conclusion can be generalized from 
the present ID case to more realistic 2D interfaces. Similarly to what was re- 
cently done for a real 2- length sliding contact (graphite flake on graphite) [2], 
our three length situation could, in principle, be accessed experimentally by 
driving, e.g., a graphite flake between two different but well defined crystalline 
surfaces. 
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a 



Fig. 1. Schematic drawing of the lubricant chain confined between the two rigid 
substrates. As shown, the sliding interface is characterized by the three inherent 
length scales: a, b and c. 
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Fig. 2. Dependence of the kinetic friction force /ki n on the chain stiffness K for 
different values of the external driving velocity. Both cases of golden and spiral 
incommensurability are shown. 
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Fig. 3. Atomic velocities versus time for (a) spiral incommensurability 
a/b = 351/265, c/a = 265/200; (b) golden incommensurability a/b = 233/144, 
c/a = 144/89. For clarity, only the velocities of 3 adjacent particles of the chain are 
shown. Panel (c) displays the corresponding Fourier transform of ii, showing the 
sharp peak exactly located at the frequency v = 1/T ~ 0.03819 that characterizes 
the time-periodic dynamics of the golden-mean case. 
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Fig. 4. Dynamics of the four sliding states marked by the black circles in Fig. 2. 
The four upper panels show the temporal behavior of the spring force, to which 
the trivial jV e ^t drift term is subtracted, with the characteristic transition from 
stick-slip motion at low K (left) to smooth sliding at high K (right). The lower 
panels display the corresponding time dependencies of the velocities of the top 
substrate and of the center of mass of the lubricant chain. 
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